REVTgX 3.0 



ANL-MSD Theory Group Preprint, 1998 



Theory of Scanning Tunneling Spectroscopy of Magnetic-Field- Induced Discrete Nodal States in a 

D-Wave Superconductor 



00 
On 
On 

o 

Q 

m 

O 
o 

Of 

3 : 
S: 

i , 

c ■ 
O ' 

o : 



> 

oo 
m 
(N 

00 
ON 
-i— > 

c3 



I 

c 

o 
o 



X 



Boldizsar Janko 

Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439 

(December 23, 1998) 

In the presence of an external magnetic field, the low lying elementary excitations of a d-wave superconductor 
have quantized energy and their momenta are locked near the node direction. It is argued that these discrete states 
can most likely be detected by a local probe, such as a scanning tunneling microscope. The low temperature 
local tunneling conductance on the Wigner-Seitz cell boundaries of the vortex lattice is predicted to show peaks 
spaced as ±-y/n, n — {0, 1, 2, ...}. The n = peak is anomalous, and it is present only if the superconducting 
order parameter changes sign at certain points on the Fermi surface. Away from the cell boundary, where the 
superfluid velocity is nonzero, each peak splits, in general, into four peaks, corresponding to the number of 
nodes in the order parameter. 



PACS numbers: 74.25.Jb, 61.16.Ch 74.60.Ec 

The nature of the quasiparticle spectrum in the Abrikosov 
vortex state of the cuprate superconductors has attracted 
strong continuing interest in recent years [jjj]. There are two 
main reasons why this spectrum is expected to be uncon- 
ventional. First, the pair size £o seems to be comparable 
to the interparticle distance l//cp. Indeed, angle-resolved 
photoemission experiments |Q| on Bi2Sr2CaCu20s+ x give 
kp ~ 0.7 A -1 , while magnetization studies [j^] indicate, that 
£o ~ 10 — 15 A. Second, there is by now substantial ex- 
perimental evidence [Q] for nodes in the cuprate supercon- 
ducting gap. In contrast to conventional superconductors, 
where the main source of low lying excitations in the mixed 
state are the bound vortex core states [Q], the above men- 
tioned reasons turn the bound core states into a delicate feature 
|^-|^] of cuprate superconductors. Experiments seem to mir- 
ror this uncertainty: while bound core states were observed in 
YBa2Cu306-<5 by far-infrared spectroscopy ^ and scanning 
tunneling spectroscopy (STS) [|To|], they seem to be absent in 
Bi 2 Sr 2 CaCu 2 8+x fy. 

The aim of the present paper is to point out, that a dis- 
crete spectrum should nevertheless be observable in the lo- 
cal density of states of the mixed state of all cuprates. It is 
somewhat unexpected, however, that such a spectrum is more 
likely to be observable close to the Wigner-Seitz cell bound- 
aries of the Abrikosov lattice, where the superfluid velocity is 
small. This is in contrast to conventional Caroli-de Gennes- 
Matricon bound states [0], which appear in the vortex core, 
in the region of singular superflow. STS is an ideal probe 
[ fl^ ] to detect these low lying discrete states, provided that the 
intervortex region is studied, and not the core region, where 
the recent STS studies [ |lO|Jl"T| ] were focused so far. The ex- 
istence of field-induced discrete quasiparticle states in a su- 
perfluid with gap nodes was sugg ested by Volovik [[TJ] for 
the A phase of superfluid 3 He, where an effective magnetic 
field is supplied by a static order parameter texture. The ana- 
logue of these nodal states in cuprate superconductors was re- 
cently proposed by Gor'kov and Schrieffer [14|, as the ap- 
propriate basis states for discussing the de Haas-van Alphen 
oscillations in the mixed state of a d-wave superconductor. 



Quite recently, Anderson [|l5|] attributed the anomalous mag- 
netothermal conductivity of Bi2Sr2CaCu20s+ x JTrj ] to the 
appearance of the magnetic field-induced discrete spectrum at 
the gap nodes. These global probes, however, cannot provide 
direct spectroscopic evidence for such a fine structure in the 
electron spectrum: spatial averaging results in averaging over 
large Doppler shifts Jl^ ] due to the supercurrents surrounding 
the vortices, which in turn act to smear out the discrete level 
spectrum [[T|]. 
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FIG. 1. (a) Schematic picture of a normal state quasiparti- 
cle momentum precessing according to <f>(t) = ui c t; (b) in the 
presence of an anisotropic order parameter, the quasiparticle mo- 
menta of the low-lying excitations are bound to the node region 

^ = (2j + l)7r/4,i = {0,l,2,3}. 

One can anticipate the existence of such a discrete spectrum 
by the following quasiclassical argument. Let us take a two 
dimensional circular Fermi surface for the quasiparticles in 
the (a,b) plane and and orient the magnetic field B along the 
c axis [see Fig. [j](a)]. The quasiparticle equation of motion in 
the normal state is 



e k 
— v k x B , 

cn t 



(1) 



where v k = h^ 1 d] l e(k) is the quasiparticle velocity. While 
the normal state quasiparticles are strongly damped, the quasi- 
particle lifetime r increases rapidly below the supercon- 
ducting transition jl9[]. If that would be the only change 
for temperatures below T c , the solution to the above equa- 
tion could be given in terms of a phase variable k = 
k (cos <f>(t), sin <j>{t)), 
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eB J FS v k 



(2) 



where oj c — (eB)/(m*c). More importantly, however, the 
onset of superconductivity generates the appearance of a gap 
tied to the entire Fermi circle [see Fig. 1(b)]. If A(</>) is 
anisotropic, such as in the case of a d-wave superconductor 
with A(<fi) = A cos(20), the order parameter will provide 
an off-diagonal confining potential in the angular coordinate 
space, for the Bogoliubov quasiparticles precessing around 
the Fermi surface. 

In order to provide a more detailed discussion of these phe- 
nomena, it is necessary to obtain the quasiparticle amplitudes 
u n ((f>), v n {4>)- Let us focus on the intermediate field regime 
B c i <C B <C B C 2 when the field in the sample is fairly uni- 
form, and address the effect of superflow later. It is most con- 
venient to use the Bogoliubov-de Gennes equations for the 
angular amplitudes of quasiparticles near the Fermi surface 
k ~ hp, in the form given by Gor'kov and Schrieffer [|l4gl: 

(E n - itkJcdj + ft)u n ((f>) + A(4>)v n (4>) = 0, 

(E n +i?kJ c d <f> - jx)v n {4>) +A(4>)u n ((f>) = 0. (3) 

where fi is defined by [i — p, = lu c Nq, (Nq is a large in- 
teger). In this approximation, the dependence of the gap 
and the amplitudes on the radial component of the momen- 
tum is neglected, and the angle <fi remains the only dynamic 
variable. The amplitudes obey periodic boundary conditions 
(u n ((f>),v n ((f>)) = (-u n (0 + 27r),w n (0 + 27r)). The normaliza- 
tion conditions are 



2tt 



\u n {4>)\ 2 + \v n {cj>f 



1. 



(4) 



As will soon become evident [cf. Eq.(JTl|)], E n ^> fi, 
and therefore ft will be neglected in the following calcula- 
tions. Using the observation of Bar-Sagi and Kuper pl"|], and 
Kosztin and coworkers [ p2[ ] for the Andreev Hamiltonian, the 
square of the Hamilton operator acting on (u„ (</>), v n (</>)) in 
(^), can be diagonalized with the help of a unitary transfor- 
mation. The eigenvalues and the eigenfunctions of eigenvalue 
problems obtained in this way 

[-h 2 cu 2 d 2 + A 2 (0) ± huj c (fyA(0))] $±, n (0) 

= A n $±, n (0) . (5) 

are simply related to the eigenvalues and eigenfunctions, 
E n , (u n {4>),v n {4>)), of the original problem (§): 



E„ 



Vn(4>) = i[$- )B (0)-**+,„(0)] 
ll 



(6) 



Furthermore, the eigenfunctions of the two branches are inter- 
related for \E„\ > : 



\E„ 



+ ,n 



\E„ 



(7) 



Here the following notation was used Q = —huj c d c f > + 
A , = h ui c dcf, + A . In the intermediate field regime 
the excitations of interest lie deep in the node region 9 = <j) — 
4>j <C 4>j, for which the d-wave potential is linear A(0) ~ 
2Ao6>. Thus, near the nodes, Eq.(Q) takes a particularly simple 
form 



*-,»(*) 



El + 2j^Ac 
Rw2 



,Ao_ 

TlLJ r 



= o, (8) 



which can be recognized as the Schrodinger equation of a 
simple harmonic oscillator, with solutions proportional to 
exp[— (A /hu; c )9 2 ]. For typical fields B = 9 T, and gap val- 
ues [g] A ~ 40meV, the ratio 7 = 2A Q /(hcj c ) ~ 80. This 
means that the significant weight of these states is strongly 
localized around the node regions. Extending the limits of 
integration to infinity in the normalization condition (Q), we 
obtain for the eigenfunctions 



*-,n(0) = C n Hn[l 1/2 0]e-l 



(9) 



where H n (x) are the Hermite polynomials, and the coefficient 



(tA) 1/2 



(5 7 /4n) + (2 7 )V2 



n > 0. 



(10) 



The angular amplitudes obeying the periodic boundary con- 
ditions can be constructed from the above solution, Eq.(^j), 
by inserting them into eqs.(||) and (0) with the substitu- 
tion <S>_,„(0) = t[mitij{^ - (2j + l)7r/4}]. The over- 
lap of the nodal states residing in two adjacent nodes is 
cxp[-(7r 2 A )/(8?kj c )] ~ 3.7 x 10~ 22 . Clearly, quantum 
oscillations arising from interference between low lying nodal 
states located at different nodes will be significant only near 
H C 2, where hu c ~ Ao. The eigenvalues are E n = AhAqTiuJc, 
which gives two branches of eigenvalues for the original equa- 
tions (if) 



E n = ±2y/ uAqTUOc 



(11) 



The spacing between these states is considerably larger than 
the magnetic energy Uui c -C 2^/A huj c <C A . For B = 9 T, 
6E n = (2 7 ) 1 / 2 ?icj c - 13 meV, while Hlu c ~ 1 meV. 

The n = state is anomalous, as it does not belong to ei- 
ther of the sets in (^). It is a so-called "zero-mode" a feature 
that is expected for the supersymmetric squared Hamiltonian 
p3|j22| ]. It is worth noting that even a strongly anisotropic 
s-wave order parameter could, in principle, give a discrete 
quasiparticle spectrum in a magnetic field. However, the con- 
ditions of normalizability required for the existence of such a 
state are naturally satisfied only for an order parameter with 
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true gap nodes [ such that A(0) ~ [d^A.^)]^—^, ((/) — <fij) 
near <f> ~ <j)j ]. 

Clearly, the experimental detection of the discrete nodal 
spectrum and especially, the observation of the anomalous 
zero mode in materials suspected of having gap nodes, would 
be quite important. Let us now find the most favorable ex- 
perimental circumstances for observing these states. As in the 
case of normal state Landau levels, the semiclassical trajec- 
tories associated with the discrete nodal states are orbits of 
radii of the order of the magnetic length d — 2eB / (he) 
fli~5| ] where B is the applied magnetic field, and he/ (2e) is the 
superconducting flux quantum. An external field of B = 9 T 
gives d ~ 150 °A. This length also gives the lattice constant 
of the Abrikosov vortex lattice: a u = d] Oa = (4/3) 1 / 4 d 
|p4|]. Since the motion occurs on the background of a strongly 
varying superfiow the spectrum is shifted [|l7|| by 5Es(r) — 
mvp ■ v s (r) typically of order SEs ~ (An£o)/(2d). This 
is not negligible when compared to the level spacing 5E n ~ 
2y/hio c A ~ 2A \/?F^/d [here^F = 2?:/k F ~ 0.23 A is 
the Fermi wavelength] . Indeed, 



SE n 
SE S 




1/2 



0.5 



(12) 



If the trajectory includes segments close to the core r ~ £0, in 
those regions the ratio drops to 5E n /5Es ~ 0.05. A global 
spectroscopic probe, such as electromagnetic absorption [§], 
will detect a spatially averaged response. Clearly, this also 
implies an average over the shift 5 Eg, which would make the 
resonant transitions between the discrete quasiparticle levels 
difficult to observe. 

The situation is quite different for a scanning tunneling mi- 
croscope, with which one could select with at least Sr ~ 1 A 
resolution the place where electrons are injected or removed 
from a superconductor. The local tunneling conductance can 
be given as 

n>0 J ~ 7r 

{\ u n{<P)\ ^ \-\VnW)\ I- ( 13 ) 



dE n 



dE n 



Here Nb — (sample area)/(27rd 2 ) is the degeneracy of 
each level, |T(0)| 2 is the square of the tunneling matrix 
element which, in principle, can be angle dependent j2^], 
E±{r, (f>, V) = E n (r, </>) ± eV and /(e) = [expQfc) + l]" 1 
is the Fermi distribution function. The local spectrum also 
depends on the angle <fr through the Doppler shift 



E n (r, 



E n ~ E s (r) cos(a + , 



(14) 



where a is the angle between the local superfluid velocity and 
the a crystal axis. For a rigid and stationary Abrikosov lat- 
tice a is fixed. Es{r) = mvFV s (r) is the maximum Doppler 
shift at particular place r. In Eq.(pj|) there is no summation 
with respect to k z , the momentum in the c direction. This two 



dimensional approximation is especially adequate for super- 
conducting Bi2Sr2CaCu2C>8+x samples p^]. 

g{V) 

2 „ , , 
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FIG. 2. The dimensionless STS tunneling spectrum g(V) 
on the Wigner-Seitz cell boundary, according to Eq.(Jl5|), for 

k B T/SE n = 1/24. 

The important observation here is that according to Eq.(Q) 
the Bogoliubov amplitudes u n ((j)),v n (4>) are highly peaked 
functions around the nodal angles. Compared to |u n (<^)| 2 and 
|i> n ((/>)| 2 , all other factors in ( |i~3| ) are weak functions of the 
angle (j> , and can be taken outside of the integral. As a conse- 
quence the integration over the angle reduces to summation 
over the nodal directions of the order parameter: For a d- 
wave superconductor, = (2j + l)7r/4, j = {0,1,2,3}. 
While the tunneling matrix element along the node directions 
\T(cj)j)\ 2 is small [p5[], it is nevertheless finite, as indicated 
by the finite zero bias conductance in tunneling experiments 
JTT|] . In the following, it is assumed that |T(0j)| 2 has the same 
value for all four nodes. The dimensionless conductance g(V) 
can be defined as 



j,n>0 I 



\ df[E-{r^V)] dflE+jr^^V)} 



d([3E n ) 



d(j3E n ) 



(15) 



In order to reveal the essential features of the above result, 
let us first assume that the tip of the scanning tunneling mi- 
croscope is positioned on the Wigner-Seitz cell of the vortex 
lattice (cf. the inset of Fig.||). For all point on the cell bound- 
ary Es(r) = 0. At low enough temperature the quasiparticle 
scattering rate [^7[|, as well as the uncertainty in SEs(r) due 
to vortex lattice fluctuations is smaller than the level spacing 



SE n > h/r qp ,5E i 



lattice 



Under these conditions the tunneling 



conductance reveals (see Fig.|2j) the discrete spectrum obtained 
in Eq.([ll]). Note that the amplitude of the zero mode is twice 
the amplitude of the finite voltage peak, since the particle and 
hole tunneling rates are the same near V = 0. This is also 
evident from Eq.( 13 ). 

As the tip is moved away from the cell boundary, towards 
one of the cores, (for example, along the dashed line of inset 
of Fig.^J), each peak splits up into four separate peaks, cor- 
responding to the four nodes of the order parameter. This is 
illustrated in the main part of Fig. [3], where, for clarity, only 
the splitting of the zero bias peak is shown. The amplitudes 
of the split peaks are four times smaller than the amplitude of 
the single peak measured at the Wigner-Seitz cell boundary. 
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FIG. 3. Main Figure: Plot of the dimensionless conductance 
g{V) [cf. Eq.(p^)] of the zero bias peak only, as the position 
is varied along the dotted line of the inset. The parameters are 
6E„ = 5meV,d = 400 A, a = tt/6, T = 6K. Inset: Schematic 
figure of a vortex lattice. The Wigner-Seitz cells, where the super- 
fluid velocity is zero, are drawn with heavy lines; the dots correspond 
to the vortex cores (not on scale with the lattice constant). The sug- 
gested path for the tip of the scanning tunneling microscope is shown 
with a heavy dotted line. 

In conclusion, this paper addressed the problem of the 
quasiparticle spectrum in a d-wave superconductor in an ex- 
ternal magnetic field. The spectrum is discrete, and the quasi- 
particle amplitudes are strongly peaked for momenta point- 
ing along the node directions. An anomalous zero mode is 
present only if the order parameter has true nodes. As it turns 
out, the level spacing is smaller than the average Doppler shift 
due to the interaction of quasiparticles with the supercurrents 
circling around the vortices. Thus, in order to detect this spec- 
trum, one needs to perform scanning tunneling spectroscopy 
on the Wigner-Seitz cell boundary of the Abrikosov vortex lat- 
tice, where the superfluid velocity is zero. The calculated local 
tunneling spectrum also indicates that slightly away from the 
Wigner-Seitz boundary, each peak splits into four subpeaks, 
as a consequence of the four different Doppler shifts arising 
for quasiparticles tunneling into the four nodes. The results 
are clearly more general than the specific case of a quasi-two 
dimensional superconductor with a d-wave order parameter, 
and can be extended in a straightforward manner to any lay- 
ered superconductor that has gap nodes. In fact, the exper- 
iment proposed here would probe intimate details of the su- 
perconducting state, such as the existence, the number and the 
position of gap nodes on the Fermi surface. 
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